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A plasma simulation technique has been developed which can be used to study the effects of 
homogeneous or inhornogeneous plasma sheath on the radiation pattern deterioration and input im- 
pedance of microwave antennas. A tank has been designed and constructed for use in the simulation 
technique which can reproduce by means of real dielectric materials the dielectric constant encount- 
ered in plasma covered antenna research. 

The radiation patterns and the input impedances of an annular slot and a thin and long rectangular 
slot have been successfully measured in the presence of a simulated lossless, homogeneous, and iso- 
tropic plasma layer of varied thickness. Comparison with the available theoretical data indicates 
a generally good agreement, although some differences exist. In the case of the radiation patterns 
these differences are attributed to the finite distance between the radiator and the receiving antenna 
on one hand and the inherent inaccuracy of the saddle point method of integration in certain regions 
on the other hand. 



1. Introduction 

Recent aerospace achievements have motivated a great number of investigations into wave 
propagation through plasma sheaths. Most of these investigations are of theoretical nature and 
they use an infinite, isotropic, homogeneous, and uniform plasma slab as a model for reentry en- 
vironment of a radiating system. In order to experimentally verify the theoretical work and to 
obtain the results for more realistic radiating system and plasma sheath configurations, a labora- 
tory simulation of a plasma sheath is very desirable. The most recent simulation attempts make 
use of artificial dielectrics, especially rodded media [Rotman, 1962; Golden and Smith, 1964; 
Golden, 1964]. These attempts are based on the fact that rodded media have refractive indices 
less than unity under certain conditions [Brown, 1953]. 

The results of these experiments generally match the theory. Certain important points, 
however, were not found to coincide with the theory. For example, the sharp peaks at the critical 
angle, sin 6 c = n, were not found as expected when an E plane slot antenna was used [Golden, 
1964]. Such discrepancies are possibly due to "higher-order reflections" from the rodded media 
which result with increasing oblique incidence [Brown, 1953]. 

It is the purpose of this paper to introduce another method of simulating the plasma sheath. 
This method will be seen to provide an exact determination of the simulated plasma parameters 
and a laboratory geometry for which it is a simple matter to accurately adjust and vary the slab 
thickness. 

The simulation method, to be described in detail in a later section, permitted the recording 
of the field patterns in the region outside of the plasma sheath. The radiators used were an 
annular slot and a rectangular slot in a ground plane. Field strength and admittance measure- 
ments were made for each antenna with various slab thicknesses. 
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2. Theoretical Background 

The problem of a long, narrow rectangular slot clad with a plasma layer attracted the attention 
of many authors. Newstein and Lurye [1956] first formulated the problem and derived the expres- 
sion for the radiation field. Omura [1962] carried out numerous calculations of the radiation 
patterns where he showed the dependence of the radiation field as affected by the thickness of 
the plasma slab and the plasma index of refraction. Tamir and Oliner [1962] interpreted the re- 
sults of the radiation field in terms of "leaky waves." Galejs [1964a, b] considered the problem 
of finding the admittance of rectangular and annular slots covered with a plasma layer. 

In the case of a long, narrow rectangular slot with a constant excitation V and covered with 
a plasma layer of thickness "a" and the equivalent index of refraction "ti" the single component 
of the far zone magnetic field is given by [Newstein and Lurye, 1956; Omura, 1962; Tamir and 
Oliner, 1962] 
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where 



G{6)= cos 0[n 2 cos cos (koavn 2 — sin 2 6)~ivn 2 — sin 2 6 sin (koavn 2 — sin 2 6)] 



(2) 



and 6 is the polar angle as shown in figure 1. 

In the case of a narrow annular slot of mean radius "6" the corresponding expression for the 
single component of the far zone magnetic field is found in the appendix to be 



H<j>o = -ia)e n 2 bVJ 1 (kob sin 0)0,6) 



e i(k r-Trl4) 



(3) 



Since the factor G(0) appears in both (1) and (3), the radiation patterns for the rectangular and 
the annular slot will be quite similar when k b < 1.5 and 6 is not too close to zero. 1 Equations 
(1) and (3) are valid for real as well as complex values of the plasma index of refraction n, i.e., loss- 
less and lossy plasmas. 



1 When k<>b < 1.5 then J x (kob) - k b. 




Figure 1. Plasma clad slot antenna geometry. 
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Figure 2. Theoretical radiation patterns of a plasma 
clad annular slot. 
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Both the numerical results of Omura [1962] and the "leaky wave" analysis of Tamir and Oliner 
[1962], as well as the simple consideration of geometric optics, predict that the radiation pattern 
will have peaks in the neighborhood of the critical angles. Moreover, as the thickness of the layer 
increases, progressively more minor peaks appear in the radiation pattern while the major lobes 
are still found in the vicinity of the critical angles [Tamir and Oliner, 1962]. This fact is also ap- 
parent in the theoretical radiation patterns of an annular slot shown in figure 2. 

The difficulty that arises in the plasma simulation attempts stems from the fact that the real 
part of the plasma dielectric constant is less than unity, namely, e p /e a < 1 where e p is the plasma 
dielectric constant and e« is the dielectric constant of free space. If, however, e p /€ a = e ps /e as < 1 
where the subscripts "ps" and "as" denote plasma simulation and air simulation, respectively, 
it is seen that a simulated plasma environment depends on the ratio of the dielectric constants 
and not on their absolute values. Thus an artificial plasma environment can be created by cover- 
ing a radiator or a scatterer under investigation with a medium having a dielectric constant less 
than that of the free space simulator. 

A major requirement of the free space simulator is that it be in the liquid form to allow move- 
ment of measuring equipment through it. Furthermore, it should have a low loss tangent to mini- 
mize signal attenuation in the medium. The plasma sheath simulator can be any foamy material 
with e r ~ 1 or simply air itself. Such a combination of dielectric materials will simulate a plasma 
with 0<e p < 1. 

With such a simulation technique, the plasma parameters can be scaled and properly defined. 
Consider the rectangular and the annular slabs as examples. If in (1) or (3) we set 



n s — n 
k s a s = k a 
k s b s = kob, 



(4) 



where the subscript "5" denotes the equivalent simulated parameters, then it will follow that the 
ratio HsIHq will be constant for all angles 6. 
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The results of (4) can also be obtained directly from Maxwell's equations [Stratton, 1941]. 
It suffices to require that the two ratios remain invariant 



C\ = jJL€ 



c 2 = fJicr- 



Z 2 



(5) 



where r is the characteristic period and / denotes length. In the present case we are dealing with 
lossless dielectrics hence only C\ need to be considered. Now the invariance of C\ requires 



/^e s /fco 2 ==//,€/ 2 to 2 , 



(6) 



which leads directly to the latter two expressions of (4). Moreover, it is to be noted that if the 
dielectric permittivities are scaled by some factor, say k, then this factor must apply to all regions 
of space. Thus, if €i refers to the plasma layer and €2 to the outer region then it follows that 

_ €1 = 6 2 

and consequently 



5i— Sn 

€2 62s 



(7) 



which is equivalent to first expression of (4). 

It should be noted that configurations involving anisotropic plasma layers such as considered 
by Wait [1961] , Hodara and Cohn [1962] , and others cannot be simulated by the present technique. 

3. Experimental System and Results 

3.1 Simulation System 

The study of a table of dielectric materials [Von Hippel, 1958] at 10 GHz has revealed that the 
requirements of the free space simulating medium, i.e., low loss-tangent, noncorrosiveness, and 
stability are satisfactorily met by Aroclor 1232. The combination of air and Aroclor 1232 (e r = 2.78 
and tan 8 = 0.008) will result in an n s = 0.60 corresponding to an electron density per cubic 
centimeter N=7.83 X 10 9 /o, where / is the actual antenna operating frequency in GHz. With the 
simulating frequency f s = 10 GHz, air layer a s = 2.9 cm, slot mean radius 6, = 0.675 cm, and the 
simulating tank containing Aroclor 1232, it follows that a/ko = 1.612 and b/ko = 0.375. As a con- 
sequence of the scaling defined in (4) and since the plasma's index of refraction is a function of 
the wave frequency, the plasma environment that can be represented by this system will depend 
on the wave frequency chosen. Table 1 shows typical physical conditions of plasma sheath antenna 
environment that can be simulated with the Aroclor 1232 as free space and air as the plasma layer 
dielectric for the case of the annular slot. 

Table 1. Plasma sheath antenna environment representable by the system 
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FIGURE 3. The simulation tank. 
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A semicylindrical tank of 22 in. inside radius and 24 in. high, as shown in figure 3, was designed 
and subsequently built to the design specifications. The wall material is plexiglass, €r = 2.59. 
The flat wall of the tank is made of 1-in. plexiglass plate and the curved wall of 0.25-in. plexiglass. 
The tank holds 78 gal. of Aroclor 1232 oil, e r = 2.78. At the operating frequency of 10 GHz, the 
tank allows a separation of the transmitting and the receiving antennas by at least 28 wavelengths. 

Since the relative dielectric constants of the tank material and the oil are so close in value, 
the interface between them has negligible effect on the wave propagation. The reflections from the 
curved plexiglass-air interface were successfully reduced to a desirable level by placing high 
performance microwave absorber against the outside wall of the tank. 

With the simulation tank located on top of a stationary bench, the ground plane could be 
situated conveniently at any distance in back of the simulation tank thus allowing the simulation of 
the plasma sheath of arbitrary thickness. The ground plane was mounted in the vertical plane 
atop a wooden cart equipped with rollers for ease of movement. The gap between the ground 
plane and the simulation tank was adjusted to the desired dimension by means of spacers and 
clamps to avoid accidental movement. The effects of the finite dimensions of both the ground 
plane and the simulated plasma layer were further alleviated by placing high performance micro- 
wave absorber around the finite boundaries of the simulated plasma layer. The absorber in this 
position acts as a matched termination of the plasma layer edges. 

The pick-up antenna consisted of a half-wave electric dipole on the end of a RG58A/U coaxial 
cable. The definition of the plane of the dipole was achieved by embedding the dipole in a plexi- 
glass strip while its feeding cable was taken out of the tank along an L -bracket also made of 
plexiglass. The L -bracket was rigidly connected to a vertical shaft, the axis of which coincided 
with the center of the aperture in the ground plane. This type of arrangement allowed approxi- 
mately ± 70° coverage in the i// direction for the largest simulated plasma layer of 2 in. in accordance 
with the slot coordinates as shown in figure 1. 

The availability of theoretical results for the radiation patterns of an annular slot and a thin, 
long rectangular slot motivated the experimental investigation of these two types of radiators. 
An annular slot was constructed using brass stock because of the ease for precision machining. 
Since an annular slot of such narrow width has a very small impedance, a matching transformer 
was necessary to match the slot impedance to the 50-fl coaxial feed line. A coaxial conical taper 
of 5.2 wavelengths was machined as an integral part of the slot piece and fulfilled the impedance- 
matching requirements. Nylon with e r = 3.02 at 10 GHz served as the dielectric of the coaxial 
conical taper. The desired slot had a mean radius ot 0.256-in. and 0.020-in. width. Figure 4 
shows the schematic of the slot and tapering section with detailed machining dimensions. The 
rectangular slot was constructed using two waveguide sections as shown in figure 5. Two brass 
pieces were machined and inserted in the guides, with the small dimension of the guide gradually 
tapering to a 0.0312-in. slot. 
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FIGURE 4. Annular slot construction. 



Figure 5. Rectangular slot construction. 



Because of the inherent difficulties encountered in the measurement of the slot antenna admit- 
tance at the aperture itself, the voltage standing wave ratio (VSWR) and the shift of the null from 
the short circuit to the loaded condition were measured at the input of the tapered transition sec- 
tion. These two parameters were measured by means of a standard slotted-section technique for 
different simulated plasma thicknesses. The knowledge of the VSWR and the null shift permits 
the calculation of the slot admittance by using the Smith Chart or the well-known transmission 
line equation. The measured value of the admittance could be transferred to any point on the 
taper had the variation of the characteristic admittance with the taper length been known which 
is not the case for the linear taper used. A short-circuit test performed on the taper indicates 
that this structure has negligible losses and, consequently, it can be considered as an ideal trans- 
former. Hence, the functional dependence of the layer thickness on the slot radiation admittance 
as measured at the input of the taper should be the same as it would be at the plane of the slot 
itself. 



3.2. Results 



The experimental antenna patterns are shown in figure 6 for the annular slot. Of importance 
to note is the occurrence of the radiation maxima near the theoretical critical angle, C = 37° and 
a low gain region near = 0°. This is entirely as predicted by the theory. Theoretical curves for 
a = 0.353 A.o and for a = 1.41 \o have been plotted and are compared with the experimental values 
in figures 6a and 6b, respectively. Figure 6c shows experimental curves for a = 0.706 X and for 
a = 1.06 \ . Very good agreement with theory is evident in figure 6a, but significant discrepancies 
are noted in figure 6b. As will be seen later, the same situation occurred in the case of the rec- 
tangular slot. 

The results of the admittance measurements for the annular slot are seen in figure 7. Refer- 
ence to theoretical curves [Galejs, 1964a] shows that the experimental curves are of the same 
general shapes as those predicted. One notable discrepancy is seen in the curve for the normal- 
ized conductance. The theory predicts approximately a 5 : 1 ratio between the maximum and mini- 
mum values of G/Yo, while the experimental curve shows approximately a 1.1:1 ratio. 

The radiation patterns of the rectangular slot are seen in figure 8. Again, note the occurrence 
of radiation maxima near the predicted critical angle, 6 C = 37°. Theoretical curves for a = 0.353 X 
and for a =1.41 k have been plotted and compared with the experimental values in figures 8a 
and 8b respectively. As observed previously, the pattern for the small sheath thickness of 0.353 X 
agrees quite well with theory, while the same significant discrepancies are noted for the large 
thickness. Figure 8c shows experimental radiation patterns for sheath thicknesses of 0.706 X 
and 1.06 K . 
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Figure 6. Experimental and theoretical radiation pat- 
terns of the plasma clad annular slot. 

In evaluating these results, it is well to keep in mind any approximations made in obtaining 
the theoretical radiation patterns. For the problem at hand, the method of saddle point inte- 
gration has been used. This method allows the representation of an integral in the form of a series 
of inverse powers of distance from the source [Brekhovskikh, I960]. It is assumed that a portion 
of the integrand, which in the present problem is the angular variation of the pattern, varies slowly. 
That is, its derivatives are sufficiently small so that the expansion may be limited to the first term. 
This approximation becomes better further from the source. 
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Figure 7. Admittance characteristics of the plasma clad 
annular slot. 
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Figure 8. Experimental and theoretical radiation pat- 
terns of the plasma clad rectangular slot. 
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FIGURE 9. Functional dependence of the minimum far- 
field distance versus sheath thickness. 
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Now, reference to figure 8 shows that for large plasma thicknesses, the angular variation 
becomes quite rapid near the critical angle. In fact, one can show that at the critical angle, the 
first derivative is proportional to the square of the layer thickness. The situation, then, is the 
following. In the case of small plasma thicknesses, the slowly varying assumption is valid for the 
radial spacing between source and pickup available to the simulation technique. As the thickness 
increases, however, the angular variation near the critical angle becomes more rapid and the cal- 
culated results become less accurate. Physical limitations prevent the increase of radial spacing. 
Thus, to obtain better agreement with the simulated patterns, it becomes necessary to include a 
second term in the saddle point approximation. At present, however, the second term has not yet 
been evaluated and, consequently, the theoretical results include the contribution of the first term 
only. 

Another major factor contributing to the discrepancies in the radiation patterns stems from the 
fact that the physical boundaries of the system are not sufficient to provide the minimum required 
far-field distance for thicknesses greater than 0.45 k . The dependence of the minimum far-field 
distance on the sheath thickness is seen in figure 9. This curve was plotted by considering the 
half power bandwidth of the radiation pattern around the critical angle. The half power band- 
widths were calculated from (2). Now from Silver [1949] we have 



Ohpbw ~ 60° I -I min, 



which yields rf/X, and subsequently the minimum far-field distance. Thus the radiation pattern 
for a = 0.353 A. is well in the radiation field, while that for a = 1.41 X is in the near field. 2 

The admittance curves for the rectangular slot are seen in figure 10. Theoretical calculations 
have been made by Galejs [1964b]. Comparison with these curves shows that the experimental 
curve for the normalized conductance agrees quite well with the available theory when the sheath 
thickness is larger than 0.5 k . The same can be said for the normalized susceptance curve. 



2 We are indebted to the reviewer for suggesting this method of determining the minimum far-held distance. 
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FIGURE 10. Admittance characteristics of the plasma 
clad rectangular slot. 



4. Conclusions 



In the present work only the simplest case of the plasma sheath environment was treated, 
i.e., the case of a lossless, isotropic, and homogeneous plasma. The simulation technique has been 
successfully tested only for this simplest case, but it is clear that it may be used for the study of 
more complex plasma sheath models. The inherent simplicity of this experimental system and 
its versatility suggest that the method is capable of simulating more realistic plasma sheath environ- 
ments which are not readily amenable to analytic solutions. 

The most serious limitation of the present technique is the fact that the method is not capable of 
simulating plasmas characterized by a negative dielectric constant. Moreover, the effective 
simulated dielectric constant of 0.36 achieved with the Aroclor 1232-air combination cannot be 
readily lowered because of the lack of suitable liquid dielectrics with e r > 2.8. However, since the 
reentry communication systems can be expected to work well above the plasma frequency, the 
above limitations may not be serious ones indeed. 



This research was sponsored by the Air Force Cambridge Research Laboratories, Office of 
Aerospace Research, USAF, under the contract AF 19 (628)-3834. 

5. Appendix 

Consider a narrow annular slot of mean radius "6" cut in an infinitely conducting plane of 
infinite extent and covered with a layer of isotropic and homogeneous plasma of thickness "a" 
and the index of refraction 4 V. A constant voltage "F" is assumed to be present across the slot. 
The horizontal plane z = a coincides with the interface of the homogeneous plasma layer and the 
surrounding medium. For convenience, reference will be given to the surrounding air medium 
(0) and plasma medium (1). 

It is clearly seen from the symmetry of the problem that there is no variation in the (^-coordinate, 

i.e., — r=0, and only the H& component of the magnetic field exists. Furthermore, it is noted that 
o<p 

the electric field at the aperture is 



E P =-V8(p-b)8(z). 
848 



(Al) 



Assuming the e~ iuit time-dependence and suppressing it throughout, the electric field components 
are given by 



1 djh 

ia)€o€ dz 



(A2a) 



and 



where H& satisfies the equation 



icoeoe p dp 



(A2b) 



_d_ 
dp 



l_d_ 



(pH*) 



d 2 H<t> 

' dz 2 



+ £ 2 //<»: 



- i2oj€ n 2 V8(p - b)8(z), 0<z<a 
,z>a. 



(A3) 



The solution to the problem can be formulated readily using Fourier-Bessel transform [Stratton, 
1941]. The single component of the magnetic field in each one of the two regions has the form 



H<t>o~- 



iVb(x)€o7l 2 



s: 



[yi sinh (yia)-fn 2 7o cosh (yia)] 



(A4a) 



//</)! 



iVbcoepn 2 
'' 2 



J -00 71 L72 sinh (y\a) + n z y cosh (yia)J 



(A4b) 



where 



:V^ 



(A4c) 



ki = 27rlki, 



(A4d) 



Changing to spherical coordinates in the configuration space and taking k = k sin /3 in the trans- 
form space, (A4a) is changed into a contour integral in the complex /3-plane. The evaluation of 
the contour integral by the saddle point method gives the far-zone field in the air region; 



H^ = — ia)e () n 2 bVJ\(k b sin 0)G(6) 



e i(kor-ir/4) 



(A5) 



where 



G(6) = cos [n 2 cos cos (k aVn 2 - sin 2 0) - iVn 2 - sin 2 sin (k a Vn 2 - sin 2 0)]" 1 . (A6) 



It is clear that (A5) and (A6) are valid for real as well as complex values of the plasma index of 
refraction, n. 

Equation (A5) gives the complete solution of the fields in the radiation zone. However, "leaky 
wave" poles, poles arising from the singularities of the contour integral and located in the improper 
Riemann sheet, may significantly affect the near and intermediate zones. These "leaky waves" 
exhibit a radial attenuation and their effect on the radiation field diminishes even when strongly 
excited. An extensive treatment of complex poles is given by Tamir and Oliner [1962]. 
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